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ABSTRACT 


The  scattering  properties  of  a  linear  array  of  parallel,  center-loaded, 
cylindrical  elements  have  been  investigated  with  the  ultimate  objective  of  ob¬ 
tain:  ig  information  about  the  character  of  the  array  irom  its  scattered  field. 
To  this  end,  a  set  of  integral  equations  for  the  currents  induced  in  the  linear 
array  illuminated  by  an  incident  plane  wave  were  derived  from  the  equations 
of  Maxwell  and  the  boundary  conditions  at  the  surface  of  the  array.  Using  a 
zero-order  approximation  to  the  form  of  the  axial  distribution  of  the  induced 
currents  in  the  array,  a  pair  of  complex  current  coefficients  were  calculated 
numerically  for  each  element  of  the  array  using  a  technique  incorporating 
ihe  set  of  integral  equations.  The  approximation  technique  gives  reasonable 
accuracy  in  the  calculation  of  the  /.'-plane,  far-zone,  scattered  field  from  the 
induced  currents,  provided  tne  electrical  half-length  of  the  elements  of  the 
array  is  less  than  5"/4  radians. 

The  scattered  field  of  an  eight -element  array  was  calculated  for  various 
conditions  of  impedance  loading  and  illumination  of  the  array. 

A  significant  result  of  this  investigation  was  the  discovery  that  the 
//-plane  scattered  field  of  a  linear  array  of  cylindrical  elements  illuminated 
by  a  plane  electromagnetic  wave  consists  of  two  factors:  a  reflection  factor 
and  an  interference  factor.  The  interference  factor  is  simply  the  complex 
array  factor  of  the  array  when  excited  with  a  uniform  amplitude  and  an 
element-to-element  phase  progression  of  2  -  ( </  a )  sin  0in(.  radians,  where 
(il  )  is  the  intereleir.ent  spacing  of  the  array  ir.  wavelengths  and  «ir  c  is  the 
angle  of  incidence  of  the  illumination.  The  reflection  factor  turns  out  lo  be 
the  //-plane  scattered  field  in  the  reflected  direction  where  the  interference 
factor  becomes  unity. 

From  the  interference  factor  we  determined  the  positions  of  the  grating 
lobes  and  the  minima  of  the  //-plane  scattering  pattern  of  the  array  for  vari¬ 
ous  plane-wave  illuminations. 

It  was  found  that  the  plane  cf  polarization  of  an  array,  the  number  of  ele¬ 
ments,  the  interelement  spacing,  and,  possibly,  the  resonant  frequency  of  the 
elements  can  be  determined  from  the  //-plane  scattering  characteristics  of 
the  passive  linear  array  of  cylindrical  elements. 
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THE  SCATTERING  OF  A  PLANE  ELECTROMAGNETIC  WAVE 
BY  A  LINEAR  ARRAY  OF  CENTER-LOADED  CYLINDERS 


INTRODUCTION 

The  majority  of  the  early  theoretical  investigations  of  electromagnetic  scattering 
involved  the  scatter ‘“g  of  a  plane  electromagnetic  wave  by  a  highly  conducting,  simple 
geometric  object,  such  as  a  sphere  or  cylinder.  Extensive  references  to  this  early  work 
can  be  fom  >i  King  and  Wu  (1).  In  recent  years  a  number  of  papers  have  appeared  in 
the  literature  which  treated  certain  simple  antennas  as  scattering  devices;  in  particular, 
the  cylindrical  antenna  has  received  considerable  attention  by  Chen  and  Liepa  (2)  and 
others.  Chen  studied  the  effect  of  central  loading  on  the  induced  current  on  a  thin  cylin¬ 
der  illuminated  by  a  plane  wave  at  normal  incidence.  By  suitable  variation  of  the  load 
impedance,  Garbacz  (3)  determined  ■  artain  antenna  parameters,  such  as  impedance  and 
power  gain,  from  measurements  of  the  scattering  cross  section  of  a  single-port  antenna 
for  fixed  frequency,  polarization,  and  antenna  orientation , 

The  literature  contains  extensive  treatment  of  the  antenna  as  a  transmitter,  that  is, 
as  a  transducer  of  current  and  voltage  at  a  terminal  pair  to  electric  and  magnetic  fields 
radiating  into  space.  The  inverse  problem  of  the  receiving  antenna,  that  is,  a  transducer 
of  incident  electric  and  magnetic  fields  into  current  and  voltage  at  a  terminal  pair,  has 
not  received  a  similar  amount  of  attention.  This  condition  oi  the  literature  stems  from 
tv.o  facts:  first,  the  mathematical  problem  involved  makes  the  analytical  study  of  the 
receiving  antenna  very  difficult  except  in  special  instances,  and,  second,  the  quantity 
almost  invariably  of  interest  in  the  receiving  rase  is  the  puwer  delivered  by  an  incident 
field  tc  the  load  connected  across  ihe  antenna  terminals.  This  quantity  can  be  obtained 
most  conveniently  from  the  transmitting  properties  oi  the  antenna  by  reciprocity  consid¬ 
erations,  without  a  detailed  knowledge  of  the  complicated  field  and  current  distributions 
over  the  surface  of  the  antenna.  If,  however,  we  are  interested  in  the  antenna  as  a  scat- 
terer,  we  are  concerned  not  only  with  the  power  delivered  to  the  load  but  also  with  the 
distribution  of  the  power  density  which  the  antenna  scatters  into  surrounding  space.  To 
find  the  scattered  electromagnetic  field,  we  need  detailed  information  about  the  currents 
induced  on  the  conducting  surfaces  of  ihe  antenna  which  radiate  the  scattered  field. 

Very  little  work  on  the  scattering  of  an  electromagnetic  wave  by  an  antenna  array 
has  appeared  in  ihe  literature.  To  improve  our  knowledge  and  understanding  of  array- 
scattering  behavior,  we  decided  to  make  a  theoretical  investigation  of  the  scattering  of  a 
plane  electromagnetic  wave  by  a  simple  linear  array  of  A  parallel,  cylindrical,  center- 
loaded  elements . 

The  approach  used  to  solve  the  steady -state  scattering  problem  is  first  to  find  an 
approximation  to  the  complex  currents  induced  in  the  elements  of  the  array  by  the  inci¬ 
dent  illumination.  Once  the  current  distribution,  in  amplitude  and  phase,  along  each  ele¬ 
ment  of  the  array  is  determined,  the  far-zone  scattered  field  radiated  by  these  currents 
is  easily  calculated  from  the  superposition  of  the  radiation  fields  of  the  individual  ele¬ 
ments  of  the  array. 
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THEORETICAL  SKETCH 
Theoretical  Model 

Consider  a  plane  electromagnetic  wave  ol  wavelength  \  with  its  electric  vector  po¬ 
larized  in  the  !  direction  incident  in  the  //  plane  plane)  of  the  linear  array  of  ,v 
center-loaded  cylinders  of  infinite  conductivity  positioned  in  space  as  shown  in  Fig.  1. 

The  electromagnetic  wave  is  incident  at  an  angle  to  the  normal  to  the  plane  of  the 
array.  The  central  impedances  zt  are  lumped,  or  without  physical  dimensions.  The 
cylindrical  elements  of  the  array  are  identical,  with  the  exception  ct  the  load  impedances, 
each  with  a  half-length  a  and  a  radius  a.  For  the  sake  of  simplicity,  we  assume  the  ele¬ 
ments  of  the  array  to  be  very  thin  cylinders  with  half-lengths  larger  than  400  times  their 
radius.  Also,  it  is  assumed  that  zva  a  is  less  than  0.0100.  In  the  theoretical  model  cho¬ 
sen,  only  the  induced  currents  that  flow  along  the  length  of  cylinders  contribute  signifi¬ 
cantly  to  the  far-zone  scattered  field;  consequently,  the  induced  currents  maintain  a  re¬ 
tarded  vector  potential  possessing  solely  a  ?  component.  From  the  symmetry  of  the 
array  and  the  uniformity  of  the  phase  and  amplitude  of  the  incident  electric  field  along 
the  entire  length  of  each  cylinder,  it  follows  that  the  induced  currents  and  their  respec¬ 
tive  vector  potentials  possess  even  symmetry  in  *. 

We  use  a  complex  notation  with  the  time -dependence  factor  rxp  <y  .  t)  of  the  incident 
wave  suppressed  in  the  treatment  of  steady-state  scattering.  The  rationalized  mks- 
couloumb  system  of  units  is  used  throughout  the  report.  All  alternating  (oscillating) 
quantities  are  complex  unless  otherwise  indicated. 


Derivation  of  the  System  of  Integral  Equations  of  the  Array  Currents 

Our  ultimate  purpose  is  to  calculate  the  scattered  field  of  the  array;  to  do  this  cal¬ 
culation  we  first  find  an  approximation  to  the  current  induced  in  each  element  of  the 
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Fig.  I  -  The  theoretical  model 
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array  by  the  incident  illumination.  The  induced  currents  must  satisfy  a  system  of  linear 
integral  equations,  which  we  shall  derive  in  the  following  pages. 

The  incident  electric  field  tangential  to  the  surface  of  the  m-th  cylinder  may  be  writ¬ 
ten  in  the  form 


'-'a  rxp(Jl  m)  . 


(1) 


where  t„  and  om  are,  respectively,  the  real  amplitude  and  phase  of  the  incident  electric 
field,  each  of  which  is  constant  ainng  the  cylinder.  The  induced  tangential  electric  field 
at  the  surface  of  the  m-th  cylinder,  maintained  by  the  induced  currents  and  charges  on  all 
of  the  cylinders  of  tne  array,  is  given  by 


f '  ,  ( ? ) 


■'l 


-  /  ■ 


(2) 


where  Am,  is  the  tangential  component  of  the  retarded  vector  potential  on  the  surface  of 
the  m-th  cylinder  due  to  all  the  array  currents  and  <t>„u  is  the  scalar  potential  on  the  sur¬ 
face  of  the  m-th  element  due  to  the  charges  on  all  of  the  cylindrical  surfaces  of  the  array. 
In  the  steady  state,  <t>riw  can  be  eliminated  by  means  of  the  Lorentz  condition  (4): 


dz 


) 


(3) 

where  p-  2-<  K  is  the  phase  constant  of  space,  is  the  radian  frequency  of  the  illumina¬ 
tion,  <•„  is  the  permittivity  of  space,  and  n0  is  the  permeability  of  space. 

Combining  Eqs.  (2)  and  (3),  we  have 


dz2 


(4) 


The  electric  field  across  the  very  small  gap  of  width  2  A  at  the  center  of  the  m-th  cyl 
inder  is  related  to  the  voltage  drop  across  the  center  load  as  follows: 


Z  /m<0>  • 

m  m 


(5) 


where  zm  is  the  complex  central  load  impedance  and  /m< o )  is  the  current  at  the  center  of 
the  m-th  element  of  the  array.  For  our  dimensionless  load,  A  must  approach  zero  in  the 
limit,  and  we  obtain  from  Eq.  (5) 


where  M ?>  is  the  Dirac  delta  function. 

With  a  perfectly  conducting  cylinder,  the  total  tangential  electric  field  vanishes  at 
the  surface  of  the  cylinder,  excluding  th»  infinitesimal  region  of  the  central  load.  This 
means 
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£  ♦  F.  .0 

ru  '•  <m 


(7) 


along  the  cylinder,  and  in  the  central  gap 


KL  =  +  rL  6<*)- 


(8) 


From  Eqs.  (1),  (4),  and  (8)  w_  obtain  ihe  following  differential  equation  for  A„u  ( z )  valid 
ovei  the  entire  length  of  the  m-th  cylinder: 


dzl 


+  /*2'W*>  :  ;'(/32,^)  jZm/m(0)M2)  -  K0  e*p  (;•('„,) 


(9) 


The  general  solution  of  the  inhomogenous  differential  equation  (9)  consists  of  the 
complementary  function  and  a  particular  integral.  The  complementary  function  is  given 
by 


[^<*>1  =  <-je><c,  cos  fiz  ♦  Cj  sin  liz)  , 


(10) 


where  (?  =  i/(M0e„)  1  2  is  the  velocity  of  propagation  of  electromagnetic  waves  in  free 
space  and  c,  and  c2  are  constants  of  integration.  The  constant  c2  is  zero,  since  ArnJ(z) 
possesses  even  symmetry  in  z  in  our  theoretical  model.  A  particular  integral  can  be 
written  as  follows: 

>]p  =  W'0  t<  «r0  //3>  <  1  -  cos  /3i)  exfdOJ  -  (Zmlm(0)  2)  sin  Hz),  for  0  <  z  <  h  ,  (lia) 

and 

[/l^fl  =  (->/<?)  [<£„//?)(  I  -  cos  0z)  txpljOj  +(Zm/m(0)/2)  sin  flz),  for  -A  <  *  <  0.  (lib) 
* p 

Equations  (11)  satisfy  the  differential  equation  (9)  for  -A  <  :  <  A.  In  order  for  Eqs.  (11) 
to  satisfy  Eq.  (9)  at  the  infinitesimal  gap  in  the  center  of  the  m-th  cylinder.  ».t  is  neces¬ 
sary  that  dA ^  dz  possess  a  jump  at  r=oofa  magnitude  j^Zjj o)  ^  which  in  turn  as¬ 
sures  that  d1Amadz'1  contains  a  delva  function  of  the  correct  magnitude  to  balance  the  im¬ 
pulsive  term  appearing  in  the  right-hand  member  of  Eq.  (9).  In  shorthand  form,  the 
general  solution  of  the  differential  equation  (9)  is  given  by: 

AmJ(  z '  -  (  -j/C)  [C,  cos  pi  ♦  (  £„//?)  (1  -  cos  >3/)  exp(/6m)  -  (Zm/m(0)  2)  sin  )|  z  |  ] .  tor  -A  <z  <  A  .  (12) 


Setting  z  -  th  in  Eq.  (12)  yields  the  following  expression  for  c, ; 

C,  =  sec  /3A  -A)  ]  -  (  l  -  cos  /3A)  e*p<  jt-m)  ♦ <Zm/m<  0)  21  sin, 3/.} 

Since  A^iz)  is  an  even  function,  the  expression  for  c,  becomes 

C,  =  sec  Ph  {yt’/t^fA)  -  <£„//?><  !  -  cos  fih)  expf/fJJ  *  lZJmiH)  2)  sin,3A}. 

From  Eqs.  (14)  and  (12)  we  obtain  the  following  equatic.i: 

-rl^fA)  =  (-/  <T)  sec  /'A  {(Zm/m(  0)  2)  sin  /i(  A  -  \z\) 

<■  (/C/l^fA)  -  ( £„  -yi )  ]  (cos  pz  -  cos  /"A)j  . 


(13) 


(14) 


iS* 


(15) 
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From  its  definition,  the  2  component  of  the  retarded  vector  potential  just  outside  the 
surface  of  the  m-th  cylinder  at  a  (lists-  ce  1  from  the  center  of  the  cylinder  is  given  ap¬ 
proximately  by  the  expression 


</&.')  Kmi{2.2‘)  dg‘ 


(16) 


where 


C*pf-//5[(  2  -  g  ')  2  t  } 

Am,  <*•*'>  z  - ' - - - r-j -  •  ,or  m  f  ' 

[(  i  -  i ' ,  2  t  (  m  -  i )  2  d’1  ] 

and 

exp  {-/;<!(  2  ;  ')  2  *  o2j  | 

A'mi  <*•*’> - - -  •  ,,,r 

i<  l  ~  2  '  )2  *  fl2] 

in  which  <1  is  the  center-to-center  spacing  of  the  cylindrical  elements  of  the  array. 

Equation  (16)  is  an  accurate  representation  of  the  vector  potential,  except  at  points 
very  near  the  ends  or  the  center  load  of  the  element.  By  substituting  Eq.  (16)  into  Eq. 
(15)  and  letting  ,  assume  all  integral  values  from  1  to  \,  we  obtain  the  following  system 
of  linear  integral  equations  foi  the  array  of  S  cer  er-loaded  cylinders: 


.v  rh 


J  /,  (  (M  ')  [A',ni  (  l,  i  ')  -  h'mii  A.  g  ')  ]  <tg '  -j  »tc  7,|(2m/m(0)  2)  »in  /J<  h  -  U  ) 

*  -  (  En  <31  »*pt  i,:m)  1  ( <<•*  t't  -  ins  /‘A)J,  m  1.2.1 . S  , 


(17) 


%  -  'e0  'o' 


120"  is  the  characteristic  resistance  of  space. 


Approximation  Technique 

The  solution  of  the  system  of  integral  equations  of  the  array  to  obtain  the  induced 
currents  in  the  cylindrical  elements  is  indeed  a  formidable  task  and  will  not  be  attempted 
here.  Instead,  we  shall  assume  the  induced  currents  can  be  represented  by  the  zero- 
order  approximation 


(,**)  lt(  (ms  rg  -  nis  .  /, )  *  sin  .*(/,  j  *  h  .  (16) 

where  <1,  and  #,  are  undetermined  complex  current  coefficients.  Note  that  Eq.  (18)  sat¬ 
isfies  the  end  conditions  that  /,  vanish  at  g  •  >> ,  independe-.l  of  finite  u,  and  .  King  (5) 
suggests  that  Eq.  (18)  is  an  adequate  approximation  for  use  in  the  calculation  of  the  far- 
zone  radiation  field  of  an  array  of  thin  cylinders  of  electrical  half-length  /.  less  the  ». 
The  scattered  or  reradiated  field  is  relatively  insensitive  to  small  errors  in  the  current 
distribution  on  the  array.  The  current  coefficients  u,  and  it,  can  be  determined  from  the 
set  of  integral  equations  and,  thus,  take  into  account  the  mutual  coupling  or  interaction 
between  the  elements  of  the  array. 
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Using  Eq.  (18)  in  Eq.  (17),  the  system  of  integral  equations  of  the  array  becomes 

h  . 

£  &  -*»<«*)]  '•.«*>  ><h  Lmi(h' 

t  ■  1  '  J 

=  -j(Zm/60)[(im(  1  -  cos  pH)  +  8m»in  /3a1  Fa(z)  +  /  (  E0/30li)  (  co.  0m  *j  n  0m)  Fc)z ) .  m  i.2.3,  ...N, 

(19) 

where 


Fc(z)  -  cm  Pz  -  cos  pH  , 
F,(z)  =  sin  /3(  A  -  1*|)  , 


«*,(<■»>  *  =  C<*'  Fei*)  , 

Lmil*)  -lmil A)  --  &<*>  f.«*>  • 

( *  )  —  I  (  COS  fll  '  -  COS  /3A)  Kmi  (2  V  ')  'll  '  , 

J-A 


Lmi(t)  =  J  sin  P(h  -  1**1)  kmf  (  2,  a  ')  dz’  . 

We  desire  to  determine  the  complex  constants  (ii  and  fB(  from  Eq.  (19)  in  such  a 
manner  as  to  cause  our  zero-order  approximation  to  the  induced  currents  to  hold  fairly 
well  along  the  entire  length  of  the  elements  of  the  array.  In  what  'allows  we  employ  a 
technique  similar  to  that  of  Chen  and  Liepa  (2)  in  the  determination  of  the  current  coeffi¬ 
cients  (l,  and  . 

For  a  given  value  of  pH ,  the  real  part  of  the  dominant  term  ( i  -  m)  of  the  sums  in  Eq. 
(19)  involving  /;£<( 2 )  and  <-£.(2)  varies  with  2  in  a  manner  like  Fc(z)  and  A’, < 2 ) ,  respec¬ 
tively,  except  near  the  center  and  the  ends  of  the*  elements.  Thus,  to  achieve  our  objec¬ 
tive,  it  seems  reasonable  to  divide  Eq.  (19)  into  the  following  two  sets  of  equations: 

N 

i  *  1  '• 

and 

N 

£  c"8  & 

I  *  I 

where  zc  and  z,  are  reference  values  in  z.  Equations  (20)  and  (21)  are  formed  by  re¬ 
spectively  equating  the  coefficients  of  Fc(z)  and  F, ( 2 )  on  opposite  sides  of  Eq.  (19).  For 
a  given  angle  of  incidence  and  value  of  /iA,  the  current  coefficients  u,  and  are  deter¬ 
mined  from  Eqs.  (20)  and  (21)  considered  as  a  single  set  of  simultaneous  linear  equations. 
For  best  accuracy  in  the  determination  of  the  constants  (1,  and  ,  we  choose  2  and  2, 
such  that  the  current  distribution  functions  Fc(  z  )  and  F,  1  2,)  remain  as  large  as  possible 
over  the  range  of  interest  in  A.  The  reference  values  zr  and  2S  are  chosen  as  follows: 


-#mj<A)]  -  S^tA)}  -  /<£„  30/3)<co»  t?m  ♦  /  m  1.2 . \ . 

(20) 

-)(Zm.60)  [oj  I  -  cob  pH )  c  Bin  3]  .  rn  1.2 .  \.  (21) 
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ic  0  for  0  <  .h  <  w  2  . 

ard 

2 c  0  aiid  2,  h  -  K  I  .  for  v  2  <  fik  <  2*  . 

When  the  elements  of  the  array  are  an  odd  multiple  of  half-wave3  in  length,  the  ex¬ 
pression  for  the  element  currents  takes  the  simple  forn. 

/,•  l /it)  l)f  co«  /it  .  (22) 

The  complex  constants  «,  are  determined  from  Eq.  (19),  which  becomes,  under  this  con¬ 
dition, 

,v 

-  V  d,Ww|(A)  (  ~jXm  60)  lim  sinifro  2)  ♦  j  ( #’0  30/>)lto*  c:m  4  j  h  in  <?m) ,  m  1,2,...,  A’,  (23) 

i  “I 

where  *  is  the  odd  number  of  half-wavelengths  in  the  element  length.  Note  that  Eqs.  (20) 
and  (21)  and  the  preceding  discussion  about  the  reference  values  of  •<  do  not  apply  to  this 
special  case. 

In  the  absence  of  experimental  or  theoretical  information  from  independent  sources 
on  the  currents  induced  in  an  array  of  linear  elements  illuminated  by  a  plane  wave,  it 
was  necessary  to  test  the  approximation  technique  used  in  this  report  to  calculate  the 
induced  currents  in  the  array.  The  array  currents  found  by  this  technique  do  not  satisfy 
the  set  of  integral  equations  of  the  array;  however,  the  approximate  currents  can  be 
checked  in  a  semiquantitative  manner.  To  achieve  this  end,  we  obtain  a  measure  of  how 
much  the  approximate  currents  fail  to  satisfy  the  set  of  integral  equations  (19)  of  the  ar¬ 
ray.  The  method  of  evaluation  of  the  approximate  tecfmique  and  the  results  obtained  for 
several  values  of  the  electrical  half-length  of  the  elements  are  presented  in  Appendix  A. 
From  the  work  of  Appendix  A,  it  was  concluded  that  the  approximate  technique  used  in 
arriving  at  the  induced  currents  gives  reasonable  accuracy  in  the  calculation  of  the  scat¬ 
tered  field  of  the  currents,  provided  the  electrical  half-length  of  the  array  elements  does 
not  exceed  .'»»  t. 

Also,  it  is  found  from  numerous  calculations  that  the  approximate  solution  to  the 
array-scattering  problem  presented  in  this  report  satisfies  the  reciprocity  theorem  of 
electromagnetic  fields  to  a  high  degree  of  accuracy. 


The  Far- Zone  Scattered  Field  of  the  Array 

We  shall  now  proceed  to  the  development  of  an  expression  for  the  far-zone  scattered 
field  of  the  array.  At  the  observation  point  P  of  Fig.  1,  which  is  at  a  distance  from  the 
array  larger  than  1000  times  the  largest  dimensior  of  the  array,  th"  2  component  of  the 
far-zone  retarded  vector  potential  of  the  current  in  the  i  -th  element  of  the  array  is  given 
approximately  by  the  expression 


go  r«p '  iPri 
ivr 


Ijiflz')  C*  p  |  »/i(  2  '  con  0  ♦  y!  n  i  n  0  n  i  n  </;)  )  dz  '  . 


(24) 


The  coordinates  r,  o,  and  locate  the  observation  point,  and  2  ’  and  y,'  ( i  -  1 )  d  are  the 

array  coordinates.  By  virtue  of  the  linearity  of  our  physical  system,  it  follows  that  the 
vector  potential  of  the  array  at  l‘  is  simply  the  vector  sum  of  vector  potentialr  at  the  ob¬ 
servation  point  of  the  individual  element  currents  of  the  array.  Therefore,  the  vector 
potential  of  the  array  at  P  is  given  by 
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-h 

exp  [;  ( <  ~  1 )  pd  s  i  n  0  s  i  n  </>]  I  /^  ( j'l  ')  exp  t  j  fit  1  to»  C)  di  '  .  (25) 

J-h 

In  the  direction  ( e ,  4>),  the  far-zone  scattered  electric  iield  of  the  array  is  related  to 
the  vector  potential  by  the  expression 

E j  -jfiCA*  »ii  o.  (26) 


M0  v  /3r) 


■i.rr 


Putting  Eq.  (25)  in  Eq.  (28),  and  using  the  fart  that  tne  element  current  is  an  even  function 
of  /V .  we  obtain  the  final  form  of  the  expression  for  the  steady-- tat e  scattered  electric 
field,  namely, 

V 

f?g  =  ;  “  sin  0  exp(-;/3r)  £ 

t  •  1 


r 


/•  (tt)  coa  (u  cos  *7)  du 


c  x  p  [;  (  i  -  1 )  fid  *  i  n  t*  a  i  n  .  (27) 


where  u  =  /3z'. 


In  the  calculation  of  the  scattered  field,  it  is  convenient  to  make  the  results  inde¬ 
pendent  of  the  distance  from  the  center  of  the  first  element  of  the  array  to  the  point  of 
observation. .  To  achieve  this  condition,  we  multiply  Eq.  (27)  by  the  reciprocal  of  the 
quantity  (60/ r)  exP(-j/3r)  to  obtain 


(Ej)  =  j  sin  6  ]T 


X  *1 


f 

LJo 


/  .  (u )  cos  (1/  co>  5)  du 


exp  [;  ( i  -  I )  ptl  sin  /sin  p\ 


(28) 


which  is  defined  as  the  normalized  scattered  electric  field. 


The  effective  scattering  coefficient  of  the  i  -th  element  of  the  array  is  defined  as 


Gi  exp ( ;  y; )  =1  lt  (u)  cos  (ucos  0)  du  ,  (29) 

o 

where  Gi  and  y(  are,  respectively,  the  real  amplitude  and  phase  of  the  scattering  coeffi¬ 
cient.  Introducing  c(  and  yit  in  Eq.  (28),  we  obtain 

f 

(Eg)  =  j  sin  6  ^  exp[/(t  -  \)  0d  sin  0  sin  4>  +  jy^  (30) 

i  =1 

It  should  be  understood  that  both  g,  and  y,  are  functions  of  fth ,  (d/h),  the  load  impedances, 
and  the  angle  of  incidence  of  the  illumination. 

Under  most  conditions  of  //-plane  scattering,  it  is  found  from  our  calculations  (Ap¬ 
pendix  B)  iisat  <.he  values  of  Gt  change  only  a  small  percentage  across  the  array,  and  the 
values  of  y,  change  approximately  linearly  across  the  array.  That  is,  y*  can  be  expressed 
in  the  form 


y*  =  y,  ' -  \)@d  sin  <t>inc  .  (31) 

where  4  nc  is  the  ingle  of  incidence  of  the  illumination  and  y ,  is  the  phase  of  the  sea  r- 
ing  coefficient,  of  the  first  element  for  non-normal  incident  illumination  or  the  average 
phase  of  the  scattering  coefficients  of  the  elements  of  the  array  in  the  case  of  normal  in¬ 
cidence.  If  in  Eq.  (30)  we  replace  6(  by  c,y,  its  average  value  taken  over  the  array,  and 


NRL  REPORT  6681 


9 


/,  by  the  expression  given  by  Eq.  (31),  Eq.  (30)  becomes  a  finite  geometric  scries  which 
ca  n  be  summed  to  yield  the  following  approximate  formula  for  the  normalized,  h  plane, 
scattered  field: 


L  »i  n  ( 


w  -*i  r.  / 

-  «•  X  l>  I  /  I  >  , 

I V  2)  J  '  L  \  ' 


-  +  ( \-  1)0  2 


)]■ 


(32) 


where  0  -  /j</u m  i  *  sin  0ini). 

^rom  extensive  calculations  we  found  that  Eq.  (32)  agrees  well  in  both  amplitude  and 
phase  with  Eq.  (30),  except  near  the  nulls  of  the  function  sin  </vv-  2)  sin  <0  2)  and  regions 
in  r/i  u,here  the  array  exhibits  resonance  phenomena.  Over  the  area  covered  by  this  in¬ 
vestigation,  we  find  that  Eq.  (32)  is  quite  helpful  in  the  exploration  of  the  general  charac¬ 
ter  of  the  scattered  field.  Equation  (32)  is  particularly  useful  and  fairly  accurate  in  the 
determination  of  the  angular  positions  of  the  major  or  grating  lobes  and  the  minima  of 
the  scattered  radiation  pattern  of  the  array. 


The  Positions  of  the  Grating  Lobes  and  the  Minima  of  the  //-Plane 
Scattering  Pattern  of  the  Array 

It  is  easily  demonstrated  from  Eq.  (32),  for  an  array  with  fixed  physical  dimensions 
and  load  impedances,  illuminated  by  a  plane  wave  of  a  given  wavelength  and  angle  of  in¬ 
cidence,  that  the  major,  or  grating,  lobes  of  the  //-plane  scattering  pattern  are  centered 
approximately  at  the  values  of  0  which  satisfy  the  relation 

0  2  -  trrrr 

or 


sin  0  -  i(2n •»  ■'•</>  "  *in  (33a) 

where  n  =  0, 1,  2, 3, ...  is  the  order  of  the  grating  lobe.  According  to  Eq.  (33a),  the  zero- 
order  grating  lobe  is  located  at  the  specular  observation  angle  ;  independent  of 

the  electrical  interelement  spacing  •</.  The  back-scattered  grating  lobes  are  centered 
at  0  =  0in<.  by  definition.  Setting  0  =  ;inr  in  Eq.  (33a),  we  obtain  the  condition  for  the  oc¬ 
currence  of  a  back-scattered  grating  lobe,  namely, 

*•»  ;in<  r  in”  ,l  ■  (33b) 

Table  1  shows  the  angles  of  incidence  where  back-scattered  grating  lobes  occur  for  vari¬ 
ous  interelement  spacings  in  wavelengths. 

Table  1 

The  Angles  of  Incidence  which  Yield  Back-Scattered  Grating  Lobes 


Anglts  of  Incidence  (degrees) 


0 

1 

2 

3 

4 

5 

6 

1/2 

0 

90 

_ 

_ 

_ 

■ 

1.0 

0 

30 

90 

- 

- 

- 

3/2 

0 

19.5 

41.8 

90 

- 

- 

1 

2.0 

0 

14.5 

30 

48.6 

90 

- 

Q  1 

5/2 

0 

11.5 

23.6 

36.9 

53.1 

90 

3.0 

0 

9.6 

15.5 

30.0 

41.8 

56.4 

90 
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Also,  from  Eq.  (32)  it  is  observed  that  the  position  of  the  minima  of  the  'lane 
scattering  pattern  are  approximately  determined  by  the  condition 

2  -  tkn  ,  *  =  I ,  Z,"3 . 

or 

■  i«  4>  =  *(2*17 /Nfid)  -  •  i »  ^i.c  .  (34) 

provided  (<p/2)is  not  zero  or  an  integral  multiple  of  v. 


THE  RESULTS  OF  A  THEORETICAL  STUDY  OF  THE 
SCATTERED  FIELD  OF  AN  ARRAY 

The  Model  Used  in  the  Investigation 

The  model  used  in  th  e  investigation  is  shown  in  Fig.  1.  The  model  involves  an  eight  - 
element  linear  array  witn  uniform  load  impedances  and  an  interelement  spacing  equal  to 
the  element  length.  The  physical  dimensions  of  the  array  are  considered  fixed,  while  the 
electrical  half-length  of  the  elements  ph  changes  linearly  with  the  frequency  of  the  inci¬ 
dent  electromagnetic  wave.  The  phase  of  the  incident  field  at  the  »-th  element  of  the  ar¬ 
ray  is  given  by 

om  -  UPd  »> »  <t>iuc  -  (35) 

The  cero-phase  reference  in  all  the  work  to  follow  is  the  phase  of  the  incident  field  at 
the  first  element  of  the  array.  In  our  scattering  calculations,  the  ratio  of  the  radius  to 
the  length  of  the  elements  Is  fixed  at  0.001,  and  we  keep  £,*  *  30  to  determine  the 
strength  of  the  incident  illumination.  The  numerical  investigation  of  the  scattered  f  field 
is  limited  in  this  report  to  tb e  tf  plane  of  the  array. 


The  Results  of  the  Investigation  of  Array  Scattering 

A  Typical  Set  of  Array  Current  Distributions — As  an  example  of  the  currents  in¬ 
duced  in  the  array,  we  present  Fig.  2,  which  shows  the  distribution  of  the  in-phase  and 
quadrature  components  of  the  axial  currents  induced  in  the  various  elements  of  the  model 
array  of  full-wave  dipoles  with  central  non  reactive  load  impedances  of  72  ohms  by  a 
plane  wave  incident  at  -20  degrees.  It  is  evident  that  both  amplitude  and  phase  distribu¬ 
tion  information  is  available  from  these  graphs. 

An  Example-  of  a  Polar  Scattering  Pattern  of  the  Array  —  Figure  3  shows  a  polar  di¬ 
agram  of  the  amplitude  of  the  scattered  electric  field  of  an  array  of  eight  half-wave 
elements  with  a  nonreactive  load  impedance  of  72  ohms  when  illuminated  by  a  plane  elec¬ 
tromagnetic  wave  incident  at  -40  degrees.  Note  the  major  scattering  lobes  in  the  re¬ 
flected  direction  of  40  degrees  and  in  the  forward  scattering  direction.  Also,  notice  the 
symmetry  of  the  pattern  about  the  line  of  the  element  centers.  In  reality,  the  scattered 
energy  in  the  forward  direction  is  slightly  less  than  in  the  reflected  direction,  because 
the  induced  surface  current  density  on  the  shadowed  half  of  the  cylinders  is  slightly  less 
than  on  the  illuminated  half  of  the  cylinders.* 


*U«ing  the  classical  formulas  (see  King  and  Wu  (1),  Eq.  13.3,  page  39)  for  the  surface 
current  density  induced  on  an  infinitely  long  cylinder  of  infinite  conductivity  by  a  plane 
electromagnetic  wave,  we  find  (for  Pa  =  0,01)  that  the  surface  current  density  at  the 
center  of  the  shadowed  half  of  the  cylinder  is  reduced  in  amplitude  by  6%  ar.d  shifted  in 
phase  by  1 1  degrees  relative  to  the  current  density  at  a  point  diametrically  opposite  in 
the  center  of  the  illuminated  half  of  the  cylinder. 
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Fig.  3  -  A  H- plane  polar  scattering  diagram  of  the  model  array 


The  Distribution  of  the  Normalized  Amplitude  of  the  Scattered  Electric  Field  in  /<  2  — 
In  Fig.  4  we  present  a  universal  graph  of  the  distribution  of  the  amplitude  of  the  normal¬ 
ized  scattered  field  of  the  eight -element  array  under  study.  The  figure  shows  the  re¬ 
flected  lobe  centered  at  ^/2  =  o  a  .&  the  first-order  grating  lobes  centered  at  4>i2  =  tv.  The 
solid  curve  represents  .'he  field  distribution  given  by  the  approximate  Eq.  (32),  while  the 
indicated  points  are  derived  from  the  more  accurate  formula  of  Eq.  (30).  The  agreement 
is  quite  good  considering  that  the  conditions  for  the  points  were  chosen  deliberately  to 
show  the  greatest  possible  disagreement  between  the  two  formulas.  Near  the  element 
resonance  at  Ph  ■  1.57  and  at  Ph  =  3,64,  where  the  reactance  of  the  load  impedance  ap¬ 
parently  resonates  the  antenna  elements,  Eq.  (32)  is  not  as  good  an  approximation  to  Eq. 
(30)  as  elsewhere  in  ph  (see  Appendix  B). 

The  effect  on  the  scattering  pattern  of  the  array  resulting  from  a  change  in  the  load 
impedance  is  portrayed  in  Fig.  5.  The  change  in  the  nonreactive  load  impedance  from  72 
ohms  to  600  ohms  has  only  a  small  effect.  The  points  calculated  from  Eq.  (30)  again  fall 
closely  on  the  solid  curve  representing  Eq.  (32). 

It  is  evident  from  our  study  that  the  normalized  solid  curve  of  Fig.  4  holds  essentially 
independent  of  the  load  impedance,  the  angle  of  incidence,  and  the  frequency  of  the  illumi¬ 
nation.  Supplemental  to  the  information  of  Figs.  4  a  1  5,  the  actual  amplitude  of  the  scat¬ 
tered  electric  field,  obtained  from  Eq.  (3f>),  in  the  reflected,  or  specular,  direction  of 
observation  (</>  =  o)  under  various  conditions  is  presented  in  Table  2. 


Fig.  4  -  A  universal  //-plane  scattering  pattern  of  the  model  array 
with  load  impedances:  7.  72  *  /'(3OO0A  -  740 /f}A)  oh™ 


! 

I 


Fig.  5  -  The  effect  of  a  change  in  load 
impedance  on  the  scattering  pattern  of 
the  model  array  of  half-wave  elements 
illuminated  by  a  plane -wave  at  normal 
incidence 


ogOL 
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Table  2 

Eg  at  the  Center  of  the  Reflected  Lobe  of 
the  Scattering  Pattern  of  the  Model  Array 


Load  Impedance 
z  (ohms) 

Eg  for  Various  V  alues  of  <£;nc 

0° 

-10“ 

-20° 

O 

O 

CO 

1 

o 

O 

N* 

1 

Bh  -  v! 4 

72  +  jO 

0.201 

0.201 

0.202 

0.203 

0.205 

150  + JO 

0.196 

0.197 

0.197 

0.198 

0.200 

300+J0 

0.182 

0.182 

0.182 

0.183 

0.184 

600  +  jO 

0.147 

0.147 

0.147 

0.148 

0.148 

72  +  1  [3CO/3A  -  (740//3  A)] 

0.116 

0.117 

0.117 

0.117 

0.118 

fih  -  77/2 

72  +  j0 

2.390 

2.36 

2.28 

2.15 

1.98 

150  +  ]0 

1.480 

1.47 

1.44 

1.39 

1.32 

300  +  j0 

0.859 

0.855 

0.845 

0.829 

0.805 

600  +  jO 

0.466 

0.465 

0.462 

0.457 

0.450 

72  +  j  [300/3A  -  (740//3A)] 

2.390 

2.36 

2.28 

2.15 

1.98 

fih  -  377/4 

72  +  jO 

0.955 

0.927 

0.799 

0.895 

0.933 

150 +)G 

0.925 

0.899 

0.767 

0.830 

0.861 

300+j0 

0.875 

0.853 

0.720 

0.736 

0.755 

600  +  i0 

0.811 

0.798 

0.675 

0.636 

0.642 

72  +  j  [3 00/5 A  -  (740//? A)] 

0.347 

0.345 

_ 

0.334 

0.358 

0.364 

0h  -  TJ 

72+10 

0.564 

0.723 

0.782 

0.796 

0.785 

150  +  10 

0.568 

0.727 

0.791 

0.806 

0.794 

300  +  10 

0.607 

0.784 

0.866 

0.887 

0.870 

600  +  10 

0.749 

1.010 

1.15 

1.19 

1.16 

72  +  j  [300/3A  -  (740//3A)] 

0.455 

0.472 

0.539 

0.563 

0.543 

,5A  =  577/4 

72  +  jO 

0.657 

0.570 

0.612 

0.607 

0.583 

150  + JO 

0.720 

0.614 

0.670 

0.665 

0.628 

300  +  jO 

0.877 

0.728 

0.808 

0.802 

0.747 

600  +  J0 

1.08 

0.875 

0.985 

0.976 

0.902 

72  +  J  [300/3/,  -  (740//3 a)] 

1.73 

1.31 

1.55 

1.53 

1.37 

The  Angular  Distribution  of  the  In-Phase  and  the  Phase -Quadrature  Components  of 
the  Scattered  Field —  In  Fig.  6  we  present  some  examples  of  the  angular  distribution  of 
the  in-phase  and  the  phase-quadrature  components  of  the  //-plane  scattered  electric  field 
in  the  half-space  on  the  illuminated  side  of  the  linear  array.  These  graphs  are  obtained 
from  Eq.  (30).  Note  the  near  symmetry  of  the  field  components  about  the  reflected  d Erec¬ 
tion  of  observation  for  normal  incidence  and  the  lark  of  r  ymmetry  exhibited  for  an  angle 
of  incidence  of  -40  degrees.  The  effect  of  changing  the  load  impedance  is  seen,  on  com¬ 
paring  Figs.  6a  and  6c,  to  result  in  merely  a  change  in  the  level  of  the  scattered  field 
without  appreciable  change  in  the  angular  distribution  of  its  in-phase  and  phase- 
quadrature  components. 
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(a)  Normal  incidence;  load  impedance  =  7Z  +  jO  ohms;  ,■/,  v 

Fig.  6  -  The  angular  distribution  of  the  in-phase  and 
the  phase -quadrature  components  of  the  //-plane 
scattered  electric  field  (Continued) 

R*  [£#  ]  -  in-phase  component  of 

,m  fo'  3  quadrature  component  of  Es. 


scattered  field  strength 
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ANGLE  OF  OBSERVATION  (DEGREES) 

(b)  Angle  of  incidence  =  -40  degrees;  load 
impedance  =  72  +  jO  ohms;  Ph  -  w/2 

Fig.  6  -  The  angular  distribution  of  the  in-phase  and 
the  phase-quadrature  components  of  the  //-plane 
scattered  electric  field  (Continued) 

fie  [  «  in-phase  component  of  Eg 

/"  [£g]  *  quadrature  component  of  E * 


ANGLE  OF  OBSERVATION  (DEGREES) 


(c)  Normal  incidence;  'r>ad  impedance  = 
600  +  jO  ohms;  h  -  2 


Fig.  6  -  The  angular  distribution  of  the  in-phase  and 
the  phase -quadrature  components  of 
scattered  electric  field 


the  //-plane 


He[f(>]  in-phase  component  of  FJ 
lm  K/]  quadrature  component  of  F" 
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The  Scattered  Field  Distribution  in  ph  —  In  this  section  we  present  graphs  which 
show  how  the  amplitude  and  phase  of  the  //-plane  scattered  electric  field  of  the  model 
array  vary  with  the  frequency  of  the  illumination  or  the  electrical  half-length  of  the  ele¬ 
ments  of  the  array.  The  direction  of  observation  of  the  far -zone  scattered  field  is  cho¬ 
sen  as  either  the  reflected  direction  or  the  back-scattered  direction.  All  the  results 
presented  in  this  section  are  calculated  from  Eq.  (30). 

In  Fig.  7a  we  show  the  characteristics  of  the  steady-state  scattered  field  observed 
in  the  reflected  direction  for  normal  incidence  and  a  nonreactive  load  impedance  of  72 
ohms.  Note  the  element  resonance  at  /3A  =  1.58  manifested  by  the  maximum  of  the 
scattered-field  strength.  In  Fig.  7b  we  present  the  Ph  characteristics  of  the  reflected 
field  for  normal  incidence  and  a  certain  type  of  reactive  load  impedance.  Again  we  ob¬ 
tain  a  maximum  of  scattered  field  at  /3A  =  1.57,  where  the  elements  resonate.  In  addition, 
a  large  resonance  peak  occurs  at  ph  =  3.64,  where  evidently  the  reactance  of  the  load 
resonates  or  tunes  out  the  reactance  of  the  impedance  of  the  elements. 

In  Figs.  7a  and  7b  and  in  the  figures  to  follow  we  find  that  the  scattered  field  be¬ 
comes  quite  weak  for  ph  smaller  than  0.40. 

In  Figs.  7c  through  71  we  show  the  field  scattered  in  the  reflected  direction  by  the 
model  array  plotted  versus  the  electrical  half-length  of  the  elements  of  the  array.  In 
Fig.  7c  the  central  load  impedance  of  the  elements  is  a  pure  resistance  of  600  ohms,  and 
the  illumination  is  incident  normal  to  the  plane  of  the  array.  The  large  load  resistance 
appears  to  have  obliterated  the  element  resonance  peak,  leaving  a  slowly  increasing  am¬ 
plitude  of  the  scattered  field  with  fih,  except  for  a  small  dip  at  ph  =  3.2.  Figure  7e  is 
similar  to  Fig.  7c,  as  might  be  expected,  since  only  the  angle  of  incidence  of  the  illumi¬ 
nation  differs  in  the  two  figures.  In  Fig.  7d  we  present  the  field  scattered  by  the  model 
array  of  eight  elements  with  nonreactive  load  impedances  of  72  ohms  when  the  illumina¬ 
tion  is  incident  at  -40  degrees.  The  scattering  characteristics  exhibited  resemble  those 
for  the  same  loading  but  with  normal  incidence,  as  shown  in  Fig.  7a.  The  element  reso¬ 
nance  appearing  at  ph  =  1.57  is  the  salient  feature  of  Fig.  7d. 


/9h 

(a)  Normal  incidence;  load  impedance  =  72  +  jO  ohms 

Fig.  7  -  The  //-plane*  scattered  electric  field 
observed  in  the  reflected  direction  as  a  function 
of  the  electrical  half-length  of  the  elements  of 
the  array  (Continued) 


RELATIVE  ELECTRIC  FIELD  STRENGTH 


NRL  REPORT  6681 


FIELD  STRENGTH 


0.4  0.8  1.2  1.6  2.0  2.4  2.8  3.2  3.6  4.0  44 

(b)  Normal  incidence;  load  impedance  = 

72  +  j(300  J*  -  740//l'A)  ohms 


(c)  Normal  incidence;  load  u. 
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(d)  Angle  of  incidence  =  -40 
degrees;  load  impedance  = 
72  f  jO  ohms 
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Fig.  7  -  The  i7-plane  scattered  electric  field 
observed  in  the  reflected  direction  as  a  function 
of  the  electrical  half-length  of  the  elements  of 
the  array  (Continued) 


RELATIVE  ELECTRIC  FIELD  STRENGTH 


O.  D.  SLEDGE 


04  0.8  12  1.6  20  24  28  12  16  40  44 

Sh 

(e)  Angle  of  incidence  =  -40 
degrees--;  load  impedance  = 

600  +  jO  ohms 
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(f)  Angle  of  incidence  =  -40  degrees;  load 
impedance  =  72  +  j  (300/3A  -  740//3A)  ohms 


Fi  7  -  The  //-plane  scattered  electric  field 
observed  in  the  reflected  direction  as  a  function 
of  the  electrical  half-length  of  the  elements  of 
the  array 
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Figure  71  shows  the  amplitude  and  phase  of  the  field  scattered  in  the  reflected  di¬ 
rection  by  the  model  ar*-ay  loaded  with  the  indicated  reactive  impedances  when  the  illu¬ 
minating  plane  wave  is  iwcident  on  the  array  at  -40  degrees.  These  characteristics 
closely  resemble  those  shown  in  Fig.  7b  over  the  interval  0.4  to  2.5  in  Ph.  Intuitively, 
the  similarity  seems  reasonable,  because  the  two  sets  of  curves  give  the  scattering 
properties  of  the  same  array  but  at  different  angles  cf  incidence  of  the  illumination. 
However,  in  the  range  extending  from  2.5  to  4.4  neither  the  amplitude  nor  the  phase 
characteristic  of  Fig.  7f  resembles  closely  the  corresponding  curve  in  Fig.  7b.  An  ex¬ 
planation  of  this  divergence  is  unavailable.  The  fairly  sharp  peak  at  ph  =  1.57  is,  as 
usual,  associated  with  the  fundamental  resonance  of  the  elements  of  the  array.  The  very 
high  and  narrow  peak  at  /3a  =  3.54  comes  about  by  virtue  of  the  reactance  of  the  load  tun¬ 
ing  out  the  reactance  of  the  impedance  of  the  elements  of  the  array  at  the  frequency  cor¬ 
responding  to  this  value  of  ph . 

We  present  in  Fig.  8  the  amplitude  and  phase  of  the  electric  field  back-scattered 
by  the  model  array  with  various  load  impedances  at  the  element  centers  when  illuminated 
by  a  plane  electromagnetic  wave  incident  at  -40  degrees.  In  all  three  sets  of  back- 
scattering  characteristics  we  have  a  large  peak  in  the  vicing,  of  ph  =  2.44,  where  a  first- 
order  grating  lobe  occurs  according  to  Eq.  (33b).  The  results  shown  in  the  figure  exhibit 
interference  b. tween  the  scattered-field  components  radiated  Ky  the  individual  elements 
of  the  array.  The  interference  manifests  itself  by  the  appearance  of  the  periodic  minima 
in  the  amplitude  of  the  field  and  the  concurrent  rather  rapid  advances  in  the  phase  of  the 
scattered  field.  The  interference  effects  naturally  are  absent  in  Fig.  7  because  in  the 
rejected  direction  the  component  scattered  fields  of  *he  elements  of  the  array  are  hearly 
in  time  phase  with  each  othe;\  In  fact,  the  general  nature  of  the  Ph  characteristics  of  the 
n- plane  scattered  field  observed  in  any  direction  can  be  explained  best  with  the  aid  of 
Eq.  (32).  From  Eq.  (32)  it  le  seen  that  the  scattered  field  for  a  given  load  impedance  and 
angle  of  incidence  consists  of  two  factors:  the  reflection  factor  VC1V  exp  [/( y,  +  n/ 2)  ]  and 
the  interference  factor  [»ia  (Aty>/2)/A>  >in  (^  2)1  «>p  \)4-  2] ,  where  0  =  Pd(»<n  4>+  sin </>inc). 

In  tbe  reflected  direction,  the  scattered  field  becomes  the  reflection  factor  only,  since  in 
this  direction  the  interference  factor  becomes  unity.  In  direction's  other  than  the  re¬ 
flected  direction,  the  effect  of  the  interference  factor  appears  along  with  the  effect  of  the 
reflection  factor.  The  magnitude  of  the  interference  factor  is  shown  in  Fig.  4  as  the  solid 
curve  plotted  versus  P  2 .  As  stated  earlier,  Eq.  (32)  does  not  give  accurate  results  near 
the  nulls  of  the  interference  factor,  as  indicated  by  the  presence  of  minima  instead  of 
nulls  in  Fig.  8.  However,  the  positions  of  the  minima  and  of  the  grating  lobes  in  Fig.  8 
almost  coincide  with  the  locations  of  the  nulls  and  grating  lobes,  respectively,  of  the  in¬ 
terference  factor.  The  amplitude  curves  of  Figs.  8a  and  8b  resemble  the  general  shape 
of  the  magnitude  of  the  interference  factor  (see  Fig.  4),  except  for  some  asymmetry  in¬ 
troduced  by  the  appropriate  reflection  factors  shown  in  Figs.  7d  and  7e.  In  Fig.  8c 
the  element  load  impedances  are  reactive,  and  the  amplitude  characteristic  does  not  re¬ 
semble  the  magnitude  of  the  interference  factor.  Figure  7f,  portraying  the  correspond¬ 
ing  reflection  factor,  accounts  for  the  distortion.  The  high  narrow  peak  at  ph  =  3.55  in 
Fig.  8c  corresponds  to  a  similar  peak  in  Fig.  7f,  where  the  reactance  of  the  load  tunes 
out  the  reactance  of  the  element  of  the  array. 

SOME  PARAMETERS  OF  A  PLANE-POL -‘RIZED  UNEAR  ARRAY  OF 
CYLINDRICAL  ELEMENTS  DEDUCED  l  ''M  ITS  SCATTERING 
PROPERTIED 

Up  to  this  point  v/e  have  made  a  study  of  the  //-plane  scattered  field  of  a  linear  array 
of  cylindrical  elements  illuminated  by  a  plane  wave  incident  in  the  //  plane  oi  the  array. 

In  other  words,  we  have  been  concerned  with  the  problem  of  determining  the  scattered 
field  given  the  character  of  the  array  and  its  illumination.  Now,  let  us  consider  the  in¬ 
verse  problem  of  the  determination  of  some  of  the  parameters  of  an  unknown  plane- 
polarized  linear  array  of  cylindrical  elements  from  the  character  of  its  //  plane  scattered 
field  and  the  known  character  of  its  incident  illumination. 
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(c)  Angle  of  incidence  =  -40  degrees;  load 
impedance  =  72  +  j  (300/3A  -  740//3A)  ohms 


Fig.  8  -  The  //-plane  back-scattered  electric 
field  as  a  function  of  the  electrical  half-length 
of  the  elements  of  the  arra 


Let  us  assume  that  the  plane  of  the  array  is  known;  then  one  can  determine  the  plane 
of  polarization  and,  hence,  the  //  plane  of  the  array  by  the  following  well-known  method. 
Illuminate  the  array  with  a  linearly  polarized  plane  wave,  incident  on  the  array  at  some 
convenient  fixed  angle  of  about  60  degrees  or  more.  Rotate  the  plane  of  polarization  of 
the  incident  wave  and  vary  the  frequency  of  the  illumination  until  a  sizable  back-scattered 
field  of  the  same  polarization  as  the  incident  field  is  obtained.  Now,  with  the  frequency 
fixed,  maximize  the  amplitude  of  the  back-scattered  field  by  further  rotation -tuning  of 
the  plane  of  polarizatic  i  of  the  incident  wave.  At  this  point,  the  plane  of  polarization  of 
the  array  is  approximately  the  same  as  that  of  the  illumination. 

Now  that  the  H  plane  of  the  array,  a  plane  orthogonal  to  the  plane  of  polarization  of 
the  array,  is  determined,  it  is  possible  to  deduce  further  information  about  the  nature  of 
the  scattering  array  from  the  knowledge  and  information  obtained  from  the  study  of  the 
scattering  properties  of  a  linear  array  presented  in  the  preceding  sections  of  this  report. 
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The  further  information  about  the  array  deducible  from  its  scattering  properties  consists 
of  the  number  of  elements  in  the  array,  the  interelement  spacing,  the  //-plane  radiation 
pattern,  when  used  as  a  uniformly  excited  transmitting  array,  and,  possibly,  the  funda¬ 
mental  resonant  frequency  of  the  elements  of  tbe  array. 

In  the  following  determinations  we  assume  an  illumination-array  setup  like  that 
shown  in  Fig.  1.  The  scattered  field  is  examined  only  in  the  H  plane  of  the  array 

As  one  continuously  increases  the  frequency  of  the  illumination  of  the  array  with  a 
fixed  angle  of  incidence  of  -40  degrees  and  &“)  for  the  load  impedance  of  each  element, 
the  amplitude  of  the  back-scattered  field  resembles  that  of  Fig.  8a.  It  follows  simply 
from  Eqs.  (33a)  and  (34)  that  the  number  of  elements  in  the  array  is  one  plus  the  number 
of  minima  lying  between  successive  grating  lobes  in  the  back- scattered  fieid.  From 
Figs.  8a  and  8b,  one  observes  that  the  scattering  array  has  eight  elements. 

To  determine  the  interelement  spacing  <J,  we  make  observations  on  the  amplitude  of 
the  back-scattered  field  at  the  f&r-zone  observation  points  P,  and  P2  of  Fig.  9.  From 
point  P,  we  illuminate  the  array  from  a  direction  <PX  degrees  from  a  reference  direction 
with  a  plane  electromagnetic  wave  and  observe  the  amplitude  of  the  back-scattered  field 
at  P,  as  the  frequency  of  the  illumination  is  scanned  through  several  grating  lobes  spaced 
uniformly  in  frequency  by  an  amount  a/,  .  At  point  P3  we  repeat  the  experiment  with  4>2 
equal  to  <pt  plus  90  degrees  and  find  the  grating  lobes  separated  by  a  frequency  increment 
A ft .  It  is  easy  to  show  from  Eq.  (33b)  in  conjunction  with  Fig.  9  that  the  element  spacing 
d  must  satisfy  the  following  transcendental  relation: 

.in*1  [C  2A/,rf]  +  •in*1  te/2A/j(/]  r  77/2.  (36) 

For  good  accuracy  the  positions  of  the  observation  points  P,  and  P2  are  chosen  so  that 
A/,  and  a/2  do  not  differ  by  more  than  50%  of  the  smaller  of  the  two  frequency  incre¬ 
ments.  The  element  spacing  can  be  found  from  Eq.  (36)  by  a  well-known  graphical 
technique. 

With  the  number  of  elements  and  their  spacing  determined,  it  is  a  simple  matter, 
from  the  theory  of  arrays,  to  find  the  //-plane  radiation  pattern  of  the  linear  array  when 
used  as  a  uniformly  excited  transmitting  ~.ray  at  a  given  frequency. 

When  one  measures  the  amplitude  of  the  field  scattered  by  the  array  in  the  reflected 
direction  as  a  function  of  the  frequency  of  the  illumination,  more  than  one  resonance 


— •!  d  !•- 


Fig.  9  -  Arrangement  for  the  determination 
of  the  element  spacing  of  the  array 
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makes  an  appearance,  because  the  load  impedance  of  the  elements  usually  exhibits  some 
reactance.  Figures  7b  and  7f  show  examples  of  this  situation.  To  determine  the  reso¬ 
nant  frequency  of  the  elements  from  these  reflection  u  characteristics  it  becomes  nec¬ 
essary  to  identify  the  elenent-resonance  peak.  Apparently  from  the  figures  cited  above, 
the  size  of  the  element-reuonance  peak  (at  pk  -  -n/i)  changes,  much  less  percentagewise 
than  the  other  resonant  peaks,  associated  with  the  reactance  of  the  load,  as  the  angle  of 
incidence  of  the  illumination  undergoes  a  change  from  0  to  -40  degrees.  With  this  crite¬ 
rion,  we  may  determine  the  resonant  frequency  of  the  elements  of  the  array,  which  is 
also  the  design  frequency  of  the  array.  Whether  this  criterion  for  ti  e  selection  of  the 
element-resonance  peak  is  reliable  for  various  kinds  of  reactive  loads  can  only  be  an¬ 
swered  by  further  investigation. 


CONCLUSIONS 

A  significant  result  of  this  investigation  is  the  discovery  that  the  //-plane  scattered 
field  of  a  linear  array  of  cylindrical  elements  illuminated  by  a  plane  electromagnetic 
wave  consists,  to  the  first  approximation,  of  two  factors:  a  reflection  factor  and  an  in¬ 
terference  factor.  The  interference  factor,  defined  as 


«■»  (A6/I/2) 

•  in  (  0/'  2 ) 


exp  [;'  (  S  -  1)^/2]  , 


where  <p  =  Pd(  •  in  <t>  +  »i"  0inc) ,  is  of  the  same  form  as  the  array  factor  of  a  uniform  lin¬ 
ear  array  oi  ,v  elements  excited  with  uniform  amplitude  and  element-to-element  phase 
progression.  In  the  reflected  direction,  where  <t>  =  -^inc  and  ^  vanishes,  the  scattered 
field  becomes  equal  to  the  reflection  factor,  since  the  interference  factor  becomes  unity. 
More  physically  speaking,  the  components  of  the  scattered  field  contributed  by  the  indi¬ 
vidual  elements  of  the  array  are  in  time  phase  in  the  reflected  direction.  The  reflection 
factor  is  calculated  from  an  approximation  technique  incorporating  the  set  of  integral 
equations  characterizing  the  array  and  the  incident  illumination.  It  follows  tliat  the  re¬ 
flection  factor  is  a  function  of  the  electrical  half-length  of  the  elements,  the  electrical 
spacing  of  the  elements,  the  load  impedance,  and  the  angle  of  incidence  of  the  illumina¬ 
tion.  With  the  aid  of  the  interference  factor  of  the  approximate  formula  tor  the  scattered 
field,  it  is  possible  to  discern  much  about  the  nature  of  the  //-plane  scattered  field,  for 
example,  the  positions  of  the  grating  lobes  and  the  minima  of  the  scattering  pattern  for 
various  plane-wave  illuminations. 


From  the  scattering  characteristics  of  a  passive  linear  array  of  cylindrical  ele¬ 
ments,  one  can  de<  ace  its  plane  of  polarization,  the  number  of  elements  in  the  array,  the 
interelement  spacing  of  the  array,  the  //-plane  radiation  pattern  of  the  array,  when  used 
as  a  uniformly  excited  transmitting  array,  and,  possibly,  the  resonant  frequency  of  the 
elements  of  the  array. 


No  doubt  the  technique  described  in  this  report  for  the  calculation  o'  the  scattered 
field  of  a  linear  array  can  be  extended  to  the  study  of  the  scattering  properties  of  planar 
arrays  of  cylindrical  elements. 
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Appendix  A 


THE  EVALUATION  OF  THE  APPROXIMATION  TECHNIQUE  OF 
CALCULATING  THE  CURRENTS  INDUCED  IN  AN  ARRAY 
BY  AN  INCIDENT  PLANE  WAVE 


To  test  the  induced  currents  calculated  using  Eqs.  (20)  and  (2i)  for  a  given  angle  of 
incidence,  array  loading,  and  value  of  0A,  we  substitute  the  current  coefficients  A,  and 
into  the  set  of  integral  equations  of  the  array,  Eq.  (19),  and  compare  the  complex  values 
of  the  right  and  left  members  of  each  equation  of  the  set  for  various  values  of  fa.  In 
Table  Al,  we  present  some  examples  of  the  comparison  of  the  right  and  left  members  of 
the  set  of  equations  (19).  The  percentage  magnitude  difference  is  defined  as  100  times 
the  ratio  of  the  absolute  value  of  the  difference  between  the  magnitudes  of  the  right  and 
left  members  to  the  average  of  the  two  magnitudes.  The  symbol  m  designates  the  equa¬ 
tion  associated  with  the  m-th  element  of  the  array. 

Graphs  of  the  percentage  magnitude  difference  and  the  absolute  value  of  the  angle 
difference,  both  averaged  over  the  array,  as  well  as  the  normalized  current  amplitude  of 
the  fourth  element  are  shown  in  Fig.  Al  plotted  versus  the  position  along  the  array  ele¬ 
ments.  It  is  assumed  for  this  discussion  that  the  norma lized-current  amplitude  distri¬ 
butions  do  not  change  appreciably  from  element  to  element  in  the  array  or  with  the  angle 
of  incidence.  For  a  given  vtlue  of  /U,  it  is  observed  from  both  the  table  and  the  graphs 
that  the  margin  by  which  the  calculated  currents  fail  to  satisfy  Eq.  (19)  is  in  general, 
greatest  for  the  larger  load  impedances  and  angles  of  incidence.  The  graphs  also  show  in 
general  that  the  margin  of  failure  to  satisfy  Eq.  (19)  becomes  large  in  the  regions  where 
the  normalized  current  arc;''*  ide  is  small  An  exception  occurs  in  the  case  of  flh  =  3.0, 
where  the  margin  of  failure  to  satisfy  becomes  fairly  large  in  the  vicinity  of  fa  =  0  for 
the  larger  load  impedances,  although  the  current  is  large  in  this  region.  Since  the  am¬ 
plitude  of  the  electromagnetic  field  radiated  by  a  given  differential  length  within  a  cylin¬ 
drical  element  var;es  linearly  with  the  amplitude  of  the  current  in  the  elementary  length, 
it  follows  that  one  should  obtain  a  more  realistic  error  criterion,  as  far  as  the  scattered 
field  is  concerned,  if  both  the  percentage  magnitude  difference  and  the  angle  difference 
are  multiplied  at  each  point  in  fa  by  the  normalized  current  amplitude  as  a  weighting 
factor.  Because  of  this  weighting  factor,  the  large  margins  of  failure  to  satisfy  Eq.  (19), 
which  occur  in  low -current  regions  of  pz,  do  not  introduce  large  errors  when  one  calcu¬ 
lates  the  scattered  field  of  an  array  with  the  numerical  approximation  technique  described 
in  this  report. 

Needless  to  say,  this  method  of  evaluation  of  our  approximation  leaves  much  to  be 
desired.  A  much  more  worthwhile  evaluation  can  be  made  by  comparing  the  calculated 
results  with  either  the  experimental  or  theoretical  results  of  other  research  workers. 
Unfortunately,  such  results  are  not  available  at  present  in  the  antenna  literature  as  far 
as  the  author  has  been  able  to  determine  from  a  limited  amount  of  literature  search. 
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Table  A1 

The  Degree  to  Which  the  Calculated  Array  Currents  /ail  to  Satisfy  the  Set  of  Integral 
Equations  of  the  Model  Array  at  Various  Points  along  the  Length  of  the  Elements 


rr 

_ 

Magnitude  Difference 
<%> 

Angle  Difference  j 
(degrees)  j 

SI 

Magnitude  Difference 
(%> 

Angle  Difference 
(degrees) 

0k  •  1.5;  0m  «  1.2;  Z  ■  72  +  JO  ohms;  j 

*  -40  degrees  | 

0k  »  3.0;  0m  a  2.5;  «  0  degrees; 

Z  a  600  +  JOohms 

1 

2 

3 

4 

5 

S 

7 

8 

1.4 

1.1 

2.2 

3.2 

3.2 

2.2 

1.4 

2.0 

5.1 

2.5 

2.5 

2.7 

3.3 

1 

2 

3 

4 

12.8 

9.6 

8.3 

7.6 

- line  of  symmetry 

2.3 

2.2 

1.5 

1.2 

3  >5 

3.0 

2.0 

0k  a  3.0;  0m  -  2.5;  *ltc 
Z-  803  +  JO  oh 

■  -40  degrees; 
ms 

0k  «  1.5;  0m  «  1.2;  ■*>!,<• 

Z  ■  72  +  ]0  oh] 

■  0  degrees; 

ms 

1 

2 

3 

4 

5 

6 

7 

8 

14.0 

16.3 

16.0 

18.0 

18.0 

16.1 

15.9 

16.6 

1.4 

2.0 

2.1 

2.1 

2.1 

2.1 

2.2 

1.0 

1 

2 

3 

4 

3.7 

6.0 

5J2 

5.5 

3.0 

2.4 

2.7 

2.6 

0k  m  1.5;  0t  m  1.2;  *,,e  •  0  degrees; 

Za  800  +  jO  ohms 

0k  a  3.0;  0t  a  0;*i.c  ■  -40degrees; 

Za  BOO  +  JOohms 

1 

2 

3 

4 

12.2 

14.3 

14.0 

13.8 

2.4 

1.8 

3.0 

2.0 

1 

2 

3 

4 

5 

6 

7 

8 

11.8 

13.2 

13.9 

13.5 

13.7 

13.5 

VI 

..2 

5.9 

7.1 

6.8 

8.7 

6.8 

8.9 

8.7 

7.5 

0k  m  1.5;  0m  m  1.2;  •Pi.c  -  -40  degrees. 

Z  *  800  ♦  JO  ohms 

1 

14.4 

8.0 

7.7 

8.5 

12.0 

13.5 

8.6 

7.7 

9.4 

5.9 

3.9 

3.2 

3.1 

5.2 

5.7 

4.4 

0k  *  3.0;  0m  a  0;  «i.e  a  o degrees; 

Z  »  72  +  J  (653)  ohms 

1 

2 

3 

4 

r  14.7 

16.7 

18.7 

18.7 

8.3 

6.9 

6.6 

6.8 

0k  m  1.5;  0m  -  1.2:  «  0  degrees; 

Z  «  72  -  J  (43.5)  ohms 

0k  a  3.0;  0m  a  0;  •  -40 degrees; 

Z  -  72  ♦  J  (653)  ohms 

1 

2 

3 

4 

3.7 

6.6 

j.2 

5.8 

3.5 

2.8 

3.1 

3.0 

1 

2 

3 

4 

5 

6 

7 

8 

12.4 

13.2 

13.0 

13.0 

13.0 

l'i.7 

13.1 

14.5 

10.4 

13.1 

14.0 

13.0 

13.3 

13.1 

13.2 

12.5 

0k  m  3.0;  0m  m  2.5;  4>i.c  »  0  degrees; 

Z  »  72  ♦  JO  ohms 

1 

2 

3 

4 

1.8 

2.9 

2.8 

2.8 

- line  of  symmetry 

10.3 

7.9* 

6.6 

6.2 

0k  «  3.0;/)*  »  2.5;  *  -40  degrees; 

Z  -  72  +  JO  ohms 

j 

9.8 

10.6 

13.6 

15.5 

12.5 

11.8 

15.3 

5.8 

18 

20.9 

22.0 

19.2 

20.0 

20.1 

20.4 

20.0 

i 

* 

I 

\ 


i 

I 
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Fig.  A 1  -  The  magnitude  difference  and  the  absolute 
value  of  the  angle  difference,  both  averaged  over  the 
array,  shown  as  functions  of  position  along  the  length 
of  the  elements  of  the  model  array  (Continued) 
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b)  jih  =  3.0;  Z  =  72  +  jO  ohms 


Fig.  A1  -  The  magnitude  difference  and  the  absolute 
value  of  the  angle  difference,  both  averaged  over  the 
array,  shown  as  functions  of  position  along  the  length 
of  the  elements  of  the  model  array  (Continued) 


AVERAGE  MAGNITUDE  DEFERENCE  (PERCENT) 
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1*0  w_ 


- MMMTUOC  OnOKNCE 


f, -40  DEGREES 

X  X 


-CURRENT 

•  ^•OOESREES) 


VOMHS- 


*®)  ^  =  3.0;  Z  =  72  +  j653  ohms 


,  DEFERENCE 


0  ‘  2‘«  1  A  '  j? 

(h)  £A  =  3.6;  Z  =  7  jO  ohms;  d>(i|t  =  0  degree 


TO  82% 


*•— 4NGLE 

\  difference 


magnitude 

difference 


(»)  P*  =  3.6;  Z  =  72  +  jO  ohms;  0j8c  =  -40degrt;es 

Fig.  A1  -  The  magnitude  difference  and  the  absolute 
V***  of  the  “8*'  difference,  both  averaged  over  the 
array,  shown  as  functions  of  position  along  the  length 
of  the  elements  of  the  model  array  (Continued)  8 
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(j)  Ph  ~  3-6;  Z  =  72  +  j874  ohms;  ^jnc  =  -40drgrees 


(k)  fih  =  3.6;  Z  =  72  +  j874  ohms;  <£in  =  0 degree 


Fig.  A1  -  The  magnitude  difference  and  the  absolute 
value  of  the  angl>"  difference,  both  averaged  over  the 
array,  shown  as  lunctions  c.f  position  along  the  length 
of  the  elements  of  the  model  array  (Continued) 
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(1)  Ph  =  2.6;  Z  -  600  +  jO  ohms;  0(nc  =  Odegree 


(m)  ph  =  3.6;  Z  -  600  +  jO  ohms;  >/jnr  =  -40degrees 

Fig.  A1  -  The  magnitude  difference  and  the  absolute 
value  of  the  angle  difference,  both  averaged  .'•■ir  the 
array,  shown  as  functions  of  position  along  the  length 
of  the  elements  of  the  .nodel  array 
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Appendix  B 

THE  AMPLITUDE  AND  PHASE  OF  THE  EFFECTIVE  SCATTERING 
COEFFICIENT  OF  THE  ELEMENTS  OF  THE  MODEL  ARRAY 


We  list  in  Table  Bl  tho  amplitude  c,  and  the  phase  y  ,  of  the  effective  scattering  co¬ 
efficient  of  the  i  -th  element  of  the  model  array,  where  i  runs  from  one  to  eight.  In  these 
tabulations  WA  means  the  deviation  from  the  average  and  DFS  means  the  deviation  from 
S,  where  s  =  pd  .i.  <£iBC  (see  £q.  (31)  of  the  text).  At  the  bottom  of  each  table  for  nor  - 
mai  incidence,  we  show  the  rats  value  of  the  percent  WA  of  the  c,  and  the  rms  value  of 
the  Dt'A  of  the  yi  both  taken  over  the  array.  In  the  case  of  incidence  at  -40  degrees,  the 
rms  value  of  the  WS  of  Ay,  ji  given. 

The  tables  bear  out  the  statement  made  in  the  text  that  the  values  of  c,  remain  ap¬ 
proximately  invariant  over  the  array  and  the  y,  approximately  linear  along  the  array. 
Tables  B1  (k)>  B1  (1),  and  Bl  (n)  exhibit  the  largest  deviation  from  this  condition  because 
they  involve  the  resonances  at/tt  -  tr/2and  near  -  3.6.  Note  that  the  element  reso¬ 
nance  ha*-,  little  effect  when  the  load  impedance  is  600  ohms  as  in  Tables  Bl  (q)  and 
Bl  (r). 


Table  Bl 

The  Amplitude  and  Phase  of  the  Effective  Scattering 
Coefficient  of  the  Elements  of  the  Model  Array 


t 

WA  Of  C, 
(%) 

y< 

(degrees) 

WA  Of  y, 
(degrees) 

Ay, 

(degrees) 

W5  Of  Ay, 

(degrees) 

(a)  7.  -  72  ♦  jO  ohms;  £A  =  1.0;  <f>iBC  *  0 degree 

1 

0.0522 

3.0 

77.0 

-0.1 

— 

— 

2 

-1.0 

76.5 

-0.6 

— 

— 

3 

-3.0 

78.0 

0.9 

— 

— 

4 

1.0 

77.1 

0.0 

— 

— 

........ 

line  of  symsM 

5 

■ 

1.0 

77.1 

1 

— 

-- 

6 

-3.0 

78.0 

0.9 

— 

— 

7 

-1.0 

76.5 

-0.6 

— 

— 

8 

0.0522 

3.0 

77.0 

-0.1 

. 

— 

rms  value 
_ 

— 

2.24 

— 

0.543 

— 

34 


(Table  Bl  continues) 
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Table  B1  (Continued) 


DFA  of  G, 

y  i 

DFA  of  y. 

DKS  Of  A y. 

(%) 

(degrees) 

(degre.s) 

(degrees) 

(degrees) 

(b)  7  =  72  +  jO  ohms;  ph  -  1.0; 

4>inc  =  -40  uegrees 

1 

0.0532 

2.1 

78.6 

- 

-76.5 

-2.8 

2 

0.0535 

■ 

2.1 

— 

-74.5 

-0.8 

3 

0.0529 

-72.4 

— 

-74.9 

-1.2 

4 

0.0523 

-147.3 

— 

-72.9 

0.8 

5 

0.0511 

-1.9 

-220.2 

— 

-74.0 

-0.3 

6 

0.0512 

-1.7 

-294.2 

— 

-71.9 

1.8 

7 

0.0503 

-3.4 

-366.1 

— 

-75.2 

-1.5 

8 

0.0520 

-0.2 

-441.3 

— 

— 

— 

rms  value 

— 

2.00 

— 

— 

— 

1.52 

(c)  z  -  300  +  jO  ohms;  Ph  =  1.0;  <t>,nc  -  0  degree 

1 

0.0417 

2.2 

56.6 

-0.5 

-- 

— 

2 

0.0404 

-1.0 

56.8 

-0.3 

— 

— 

3 

0.0403 

-1.2 

58.1 

1.0 

— 

4 

0.0408 

0.0 

57.1 

0.0 

— 

— 

line  of  symmi 

rms  value 

— 

1.35 

- 

0.579 

1 

— 

(d)  7.  -  300  +  jO  ohms;  P a  =  1.0;  4>inc  =  -40 degrees 

1 

0.0427 

3.4 

57.6 

— 

-76.0 

-2.3 

2 

0.0422 

2.2 

-18.4 

— 

-74.1 

-0.4 

3 

0.0417 

1.5 

-92.5 

— 

-74.4 

-0.7 

4 

0.0411 

-166.9 

— 

-72.7 

1.0 

5 

0.0406 

-1.2 

-239.6 

— 

-73.9 

-0.2 

6 

0.0406 

-1.2 

-313.5 

-72.1 

1.6 

7 

0.0404 

-..7 

-385.6 

— 

T5.5 

-1.8 

8 

0.0411 

■ 

-461.1 

— 

— 

— 

rms  value 

— 

1.75 

— 

- 

- 

1.35 

(Table  B1  continues) 
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Table  B1  (Continued) 


DFA  Ol  G( 

DFA  of  y, 

Ayi 

DKS  of 

i 

(%) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(e)  z  =  600  +  jO  ohms;  ph  =  1.0;  =0  degree 

1 

0.0294 

1.4 

45.0 

-0.6 

— 

— 

2 

0.0287 

-1.0 

45.4 

-0.2 

— 

— 

3 

0.0288 

-0.7 

46.3 

0.7 

— 

— 

4 

0.0290 

0.0 

45.0 

0.0 

— 

— 

ine  of  symm< 

_ _ 

211  J - 

rms  value 

— 

0.929 

— 

u.472 

— 

- 

(f)  z 

=  600  +  jO  ohms;  Ph  =  1.0 

;  ^inc  =  “40  degrees 

1 

0.0299 

2.8 

45.6 

— 

-75.3 

-1.6 

2 

0.0296 

1.7 

-29.7 

— 

-73.9 

1 

3 

0.0293 

0.7 

-103.6 

— 

-74.0 

4 

0.0289 

-0.7 

-177.6 

— 

-72.9 

0.8 

5 

0.0287 

-1.4 

-250.5 

— 

-73.9 

-0.2 

6 

0.0287 

-1.4 

-324.4 

— 

-72.5 

1.2 

7 

0.0287 

-1.4 

-396.9 

— 

-75.1 

-1.4 

8 

0.0290 

-0.3 

-472.0 

— 

— 

— 

rms  value 

— 

1.49 

— 

- 

— 

0.845 

(g)  z  =  72  +  jO  ohms;  P h  =  3.6;  £inc  =  Odegree 

1 

1.4 

-63.1 

-4.4 

- 

— 

2 

-1.7 

-59.5 

-0.8 

- 

— 

3 

-0.6 

-56.8 

1.9 

— 

— 

1.0 

-55.5 

3.2 

- 

— 

Line  of  symm 

etry  -  - - 

rms  value 

— 

1.25 

— 

2.91 

- 

(Table  B1  continues) 
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Table  B1  (Continued) 


i 

Df-'A  of  G{ 
(%) 

(degrees) 

Dr'A  of  y. 

* 

(degrees) 

Ay,- 

(degrees) 

CFS  of  A >■, 

(degrees) 

(h)2 

=  72  +  jO  ohms;  Ph  =  3.6; 

0inc  =  -40  degrees 

1 

0.0817 

-3.7 

-71.9 

— 

-263.6 

-2.9 

2 

0.0873 

3.0 

-335.5 

— 

-267.6 

3.7 

3 

0.0867 

2.2 

-603.1 

-264.6 

-0.4 

4 

0.0853 

0.6 

-867.7 

— 

-265.8 

1.0 

5 

0.0851 

0.4 

-1133.5 

— 

-264.5 

6 

0.0824 

-2.8 

-1398.0 

— 

-262.7 

-3.8 

7 

0.0826 

-2.6 

-1660.7 

— 

-262.9 

-4.0 

8 

0.0871 

2.7 

-1923.6 

— 

— 

— 

rms  value 

— 

2.50 

- 

- 

— 

1.74 

(i)  Z  =  600  +  jO  ohms;  Ph  =  3.6;  £inc  =  0 degree 

1 

0.128 

3.4 

-38.1 

-6.6 

— 

— 

2 

0.121 

-2.3 

-32.9 

-1.4 

— 

- 

3 

0.122 

-1.5 

-28.7 

2.8 

— 

- 

4 

0.124 

0.2 

-26.4 

M 

- 

— 

...... _ 

ine  of  symm< 

rms  value 

— 

2.18 

— 

4.45 

— 

— 

0)  z 

=  600  +  jO  ohms;  Ph  =  3.6;  <*inc  =  -40  degrees 

1 

0.132 

-51.5 

_ 

-262.3 

1.6 

2 

0.144 

-313.8 

— 

-268.9 

-2.4 

3 

0.145 

-582.7 

— 

-264.8 

0.6 

4 

0.141 

1.4 

-847.5 

— 

-266.2 

-0.6 

5 

0.141 

1.4 

-1113.7 

— 

-265.0 

0.7 

6 

0.135 

-2.9 

-1378.7 

— 

-261.4 

2.5 

7 

0.132 

-5.0 

-1640.1 

— 

-261.2 

2.3 

B 

0.142 

2.2 

-1901.3 

_ 

— 

- 

rms  value 

— 

3.51 

. 

— 

1 

1 

- 

2.77 

(Table  B1  continues 
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Table  B1  (Continued) 


i 

G , 

Dfr'A  Of  C, 
(%) 

(degrees) 

DFA  Of  >, 

(degrees) 

A>. 

^degrees) 

DFi  Of  A  y. 

(degrees) 

(k)  z  =  72  +  j  [300/3*  .  (740//3A)]ohms;  M  =  3,6;  <pjn 

t  =  0  degree 

line  of  symmetry 


rms  v^lue 


(1)  z  =  72  +  j [<5uO/3A  -  (740//3A)]onms;  fth  -  3.6;  tf>inc  =  -40degrees 


-258.4 

-273.2 

-264.9 

-267.1 

-266.4 

-258.4 

-254.8 


rms  value 


(m)  z  =  72  +  jO  ohms;  (ih  -  v. 

2;  4>itiC  =  0  degree 

1 

0.278 

-7.0 

-7.1 

-4.3 

- 

— 

2 

0.313 

4.7 

-0.5 

2.3 

- 

- 

3 

0.300 

0.3 

•2.2 

0.6 

- 

— 

4 

0.305 

2.0 

-1.6 

1.2 

— 

— 

1 

rms  value 


(Table  B1  continues) 
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Table  B1  (Continued) 


DFA  Of 

>, 

DFA  Of 

A>t- 

DFS  Of  A>, 

(%) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(n)  Z  -  72  +  jO  ohms;  flh  =  w/ 2 

;  <fiinc  =  -40 degrees 

1 

0.289 

16.5 

-21.2 

— 

-110.7 

5.0 

2 

0.233 

■|||| 

-131.9 

— 

-109.5 

6.2 

3 

0.240 

-241.4 

— 

-115.5 

0.2 

4 

0.258 

-356.9 

- 

-118.8 

-3.1 

5 

0.265 

6.8 

-475.7 

— 

-118.7 

6 

0.255 

2.8 

-594.4 

— 

-115.3 

7 

0.233 

-6.0 

-709.7 

— 

-104.8 

10.9 

8 

0.213 

-14.1 

-814.5 

— 

— 

— 

rms  value 

- 

8.83 

— 

— 

5.24 

(o)  Z  -  300  +  jO  ohms;  flh  -  ?i 

/2;  «  0  degree 

1 

0.105 

-2.0 

-3.0 

-1.95 

— 

— 

2 

0.109 

2.0 

0.3 

1.35 

— 

— 

3 

0.107 

0.0 

-1.0 

0.05 

— 

— 

4 

0.108 

1.0 

-0.5 

0.55 

— 

— 

rms  value 

— 

1.50 

— 

1.22  t 

_ 1 

_ 

(p)  Z 

=  300  +  jO  ohms;  ph  -  <pinc  =  -40  degrees 

1 

0.109 

8.4 

-8.6 

— 

-114.7 

1.0 

2 

0.0986 

-1.6 

-123.3 

-112.5 

3.2 

3 

0.0974 

-3.6 

-235.8 

-114.4 

1.3 

4 

0.101 

0.4 

-350.2 

-117.0 

-1.3 

5 

0.104 

3.4 

-467.2 

-118.0 

-2.3 

6 

0.104 

3.4 

-585.2 

— 

-116.5 

-0.8 

7 

0.0985 

-2.6 

-701.7 

“ 

-111.4 

4.3 

8 

0.0932 

-7.6 

-813.1 

— 

— 

— 

rms  value 

- 

4.66 

— 

— 

— 

2.36 

(Table  B1  continues) 
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Table  B1  (Continued) 


DKA  of  (,\ 
(%) 

7. 

(degrees) 

DFA  of  y, 
(degrees) 

Ay, 

(degrees) 

DFS  of  Ay, 
(degrees) 


rms  value 


(q)  l  =  600  +  jO  ohms;  flh  -  2;  <*>inr  =  0 degree 


line  of  symmetry 


(r)  z  =  600  +  jO  ohms;  /3A  =  77/2;  4>inr  =  -40  degrees 


-115.4 

-113.9 

-114.6 

-116.3 

-117.2 

-116.4 

-113.4 


0.0589 

4.65 

-4.8 

0.0558 

-0.85 

-120.2 

0.0551 

-2.1 

-234.1 

0.0561 

0.32 

-348.7 

0.0574 

1.99 

-465.0 

0.0574 

1.99 

-582.2 

0.0558 

-0.85 

-698.6 

0.0538 

#4.41 

-812.0 

rms  value 
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S  ASSTRAC  T 

The  scattering  properties  of  a  linear  array  of  parallel,  center-loaded,  cylindri¬ 
cal  elements  have  been  investigated  with  the  ultimate  objective  of  obtaining  informa¬ 
tion  about  the  character  of  the  array  from  its  scattered  field.  To  this  end,  a  set  of 
integral  equations  for  the  currents  induced  in  the  linear  array  illuminated  by  an  in¬ 
cident  plane  wave  were  derived  from  the  equations  of  Maxwell  and  the  boundary 
conditions  at  the  surface  of  the  array.  Using  a  zero-order  approximation  to  the 
form  of  the  axial  distribution  of  the  induced  currents  in  the  array,  a  pair  of  complex 
current  coefficients  were  calculated  numerically  for  each  element  of  the  array  using 
a  technique  incorporating  the  set  of  integral  equations.  The  approximation  technique 
gives  reasonable  accuracy  in  the  calculation  of  the  //-plane,  far-zone,  scattered  field 
from  the  induced  currents,  provided  the  electrical  half-length  of  the  elements  of  the 
array  is  less  than  5-n/A  radians. 

The  scattered  field  of  an  eight-element  array  was  calculated  for  various  condi¬ 
tions  of  impedance  loading  and  illumination  of  the  array. 

A  significant  result  of  this  investigation  was  the  discovery  that  the  //-plane 
scattered  field  of  a  linear  array  of  cylindrical  elements  illuminated  by  a  plane 
electromagnetic  wave  consists  of  two  factors:  a  reflection  factor  and  an 
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interference  factor.  The  interference  factor  is  simply  the  complex  array  factor  of 
the  array  when  excited  with  a  uniform  amplitude  and  an  element-to-element  phase 
progression  of  2v (d  k)  sin  </>inc  radians,  where  (d/k)  is  the  interelement  spacing  of 
the  array  in  wavelengths  and  4>inc  is  the  angle  of  incidence  of  the  illumination.  The 
reflection  factor  turns  out  to  be  the  //-plane  scattered  field  in  the  reflected  direc¬ 
tion  where  the  interference  factor  becomes  unity. 

From  the  interference  factor  we  determined  the  positions  of  the  grating  lobes 
and  the  minima  of  the  //-plane  scattering  pattern  of  the  array  for  various  plane - 
wave  illuminations. 

It  was  found  that  the  plane  of  polarization  of  an  array,  the  number  of  elements, 
the  interelement  spacing,  and,  possibly,  the  resonant  frequency  of  the  elements  can 
be  determined  from  the  //-plane  scattering  characteristics  of  the  passive  linear 
array  of  cylindrical  elements. 
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